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ABSTRACT 



Vortex-type motions have been measured by tracking bright points in high-resolution observations of the solar photosphere. These 
small-scale motions are thought to be determinant in the evolution of magnetic footpoints and their interaction with plasma and 
therefore likely to play a role in heating the upper solar atmosphere by twisting magnetic flux tubes. We report the observation of 
magnetic concentrations being dragged towards the center of a convective vortex motion in the solar photosphere from high-resolution 
ground-based and space-borne data. We describe this event by analyzing a series of images at different solar atmospheric layers. By 
computing horizontal proper motions, we detect a vortex whose center appears to be the draining point for the magnetic concentrations 
detected in magnetograms and well-correlated with the locations of bright points seen in G-band and CN images. 
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1. Introduction 

Convection is an important mechanism of energy exchange re- 
sponsible for the formation of the extensively studied solar gran- 
, ulation. The interaction between the convective plasma flows 
' and solar magnetic fields causes the appearance of solar struc- 
tures on many spatial scales, sunspots and pores being the most 
^ conspicuous ones. Magnetism is often understood to be the 
i/^ ■ main process controlling the whole Sun's configuration and be- 
00 ! havior, even in the so-called quiet regions where small mag- 
^— H netic elements such as bright points (hereafter BPs) are de- 
tected mainly along dark lanes in-between granules. Small-scale 
Tu-jl (<0.5 Mm) convectivel y driven vortex-type motions were dis- 
■ co vered by iBonet et al.l ((2 008). These motions were suggested 
I by ISturrock & Uchidal ([^81) to account for about half of the 
photospheric kinetic energy, but this could not be confirmed 
due to the limitations of early observations. On larger spatial 
, , . scales, i.e., g ranular sizes, vortical motions were observed by 
^ ^ iBrandt et al.r(ll988h . lTarbeU et al.1 ([T99Th . and lTitle et"aD (Il992h . 
?-H Signatures of photospheric vortex flows have also been ob- 
. 5^ served on even la rger scales such as supergranular junctions 
(lAttie, Innes & Pottsl i2009). Vortical m otions in the solar pho- 
tosphere were studied by simulations ( Zirker |[l993l and refer- 
ences herein) and are understood to play an important role in the 
evolution of magnetic footpoints and heating in the upper so- 
lar atmosphere (Brandt et al. 1988; Wang et al. 1995). Here we 
describe a vortex motion in the solar photosphere that persists 
for about 20 minutes. It seems to affect the small-scale magnetic 
concentrations, which appear to be dragged towards the vortex 
center. 



2. Observations and data processing 

An international campaign at the Canary Islands observato- 
ries was performed during September- October 2007 and corre- 
sponds to an important collaboration among several European 
and Japanese institutions. Coordinated observations involving 
all the ground-based fa cilities at Canary Is lands and the space 
solar telescope Hinode (iKosugi et al. I l2007h were performed as 
part of the Hinode Operation Program 14. In the present work, 
we use data recorded during a particular observing run on 
September 29, 2007, whose target was a quiet sun region north 
of the active region NOAA 10971, close to the solar disk center 
Ou=0.99). 



2.1. Ground-based observations 



We used the Swedish 1-m Solar Telescope (SST jScharmer et aP 
.2003.) to acquire images in G-band (/1430.56 nm) with an 
eflTective field-of-view (hereafter FOV) of 68^54x68^54 and 
a sampling of 0.034"/px. Image post-processing included the 
Multi-Object Miilti- Frame Blind-Deconvolution (MOM 
I Van Noort. Rouppe van der Voort & Lofdahll l2005h restoration 
technique to correct for atmospheric and instrumental aberra- 
tions. The final products were two time series of images (si: 
08:47-09:07 UT and s2: 09:14-09:46 UT) with a cadence of 
15 s. Additional data post-processing steps were: compensation 
for diurnal field rotation, rigid alignment of the images, correc- 
tion for disto rtion, and s ubsonic filtering to eliminate residual 
jittering (Title et al. Ill 986h . 



2.2. Satellite data 



Send offprint requests to: L. Balmaceda, e-mail: 
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The Solar Optical Telescope (SOT, iTsuneta et aP 12008") on- 
board Hinode acquired filtegrams in CN (/1388.3 FWHM 0.8 
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Fig. 1. Quasi- simultaneous observations showing the evolution of a quiet sun region on 29 September 2007. From top to bottom: Mg 
I magnetogram, G-band, CN, Ca II H. Every panel displays the average image over 4-min intervals and the time stamps correspond to 
the initial time for each interval, respectively. BPs are detected in all filtergrams and are co-spatial with the magnetic concentrations 
observed in magnetograms. Spatial scales in both axis are in arcsec and correspond to solar coordinates (FOV ~6''x5'0- The black 
box at 09:04 UT shows the FOV covered by the SP in Fig.[2l Encircled in white in G-band frames at 08:48, 08:52, and 08:56 UT is 
the location of the detected vortex (Sect. 13.31) for reference. 



nm) and the core of Ca II H line (/1396.8 ± 0.3 nm) using the 
Broadband Filter Imager (BFI) with a cadence of ~35 s. The ob- 
served FOV covered 19^.^8 x 74709. Although Hinode's FOV 
is smaller than that of SST, they still overlap. Magnetograms 
in the Mg I line (/1517.3 nm) were also obtained at a cadence 
of -20 s with the Narrowb and Filter Im ager (NFI). Data from 
the spectropolarimeter (SP, llchimoto et al. 2008) were also ac- 
quired. This dataset consisted of Stokes parameters I, Q, U, and 
V measured along a slit of 256 pixels in raster scan mode of 
18 scans from 08:20-09:44, with a sampling of 0.15'7px and 
FOV of 2766 X 4a.'57. The operation mode was set to dy- 
namic mode, i.e., the exposure time per slit position was 1.6 s 
to study the evolution of highly dynamic events. The noise level 
was 1.6 X 10"^/c in Stokes I and 1.8 x 10"^/^ Stokes Q and 
U. The observed lines were the two Fel lines, /1630.15 nm and 
/1630.25 nm. The SOT images were corrected for dark current, 
flat-field, and cosmic rays with standard procedures. A subsonic 
filter was applied to remove high frequency oscillations in fil- 
tergrams. SP data were inv erted using the full atmosphere inver- 
sion code LILIA (|Socas N avarro 2001), based on the SIR code 
(iRuiz Cobo & Del Toro Iniestaiil992i) . 

3. Data analysis and results 

3.1. Multiwavelength description of the event 

Quasi- simultaneous data taken by Hinode and SST were used. 
The maximum time-lag between CN and Ca images was ~2 s 
and ~12 s, respectively, between filtergrams and Mgl magne- 
tograms. Images were rotated, co-aligned, and trimmed to main- 
tain the same FOV at all wavelengths. Figure [T] shows from bot- 
tom to top, time-averaged Ca II H, CN, and G-band filtergrams 
and Mg I magnetograms from 08:40 to 09:20 UT. Each frame 
corresponds to an average over 4 min, obtained from ~7 images 
in the case of filtergrams and ~12 Mgl magnetograms. The av- 
eraging of the consecutive images allows us to easily identify 
the structures present in successive images. We superimpose on 
the corresponding filtergrams magnetogram contours, which are 
useful for visualizing the co-spatiality of the features observed 
at different wavelengths. The grayscale bars indicate the inten- 
sities normalized to the mean intensity of the entire data set for 
the respective wavelengths. The LOS magnetograms shown here 



consist of masks obtained after applying a threshold of 3cr to re- 
move the noise following iKrivova & Solankj (|2004|) . The good 
agreement between the magnetic features detected in these im- 
ages and the bright elements present in the filtergrams gives us 
confidence to justify its use. On first inspection, the magnetic 
feature appears to undergo an apparent rotation over the 40 min- 
utes elapsed. We can identify a magnetic structure with two main 
lobes (denoted by LI and L2 in the first frame, top row of Fig. [B 
one of which rotates around the other as indicated by the plus 
black and white plus-signs marking the 4-min averaged location 
of the centroid for LI and L2, respectivel}Q. This movement is 
also detected at the various wavelengths, i.e., diff'erent heights 
in the solar atmosphere. Both G-band and CN data provide in- 
formation about the photospheric level, while Ca II H images 
and Mg I magnetograms retrieve information from the low chro- 
mosphere. During the interval between 09:00 and 09:08, the L2 
lobe seems to dissolve. The magnetic flux density in this region 
appears to decrease, but increases again when the whole struc- 
ture has rotated and L2 approaches the lower part of the FOV. An 
increase in intensity is observed in Ca images and is accompa- 
nied by the appearance of BPs in G-band, while in CN the signal 
is more difl'use. As a consequence of the rotation, small-scale 
processes of fragmentation and coalescence of BPs clearly dis- 
cernible in G-band filtergrams take place along the intergranular 
lane. The BPs are also detectable in CN images, although their 
appearance is not so well defined as in the SST data, especially 
towards the end of the sequence. The good agreement between 
the small elements detected in Ca II H and the other wavelengths 
is also remarkable. A forthcoming work (Vargas Dominguez et 
al., 2010; hereafter Paper I) will focus on the evolution of BPs in 
the sequence. 

3.2. Characterization of the magnetic configuration of the 
region 

Figure [2] displays the evolution of the magnetic structure in the 
sequence of single maps of magnetic field strength \B\ retrieved 
by LILIA (first row) and LOS-velocity, Vlos , obtained by esti- 



^ To identify the lobes, we use the line orthogonal to the linear fit of 
the pixels belonging to the structure and that passes through the center- 
of-mass of the whole data cube. 
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Fig. 2. From top to bottom: magnetic field strength and LOS velocity single maps obtained by inverting SP data. The FOV covered is 
shown in the left lower panel. Positive/negative Vlos values correspond to upward/downward velocity directions. Encircled region 
in Vlos map at 08:44 UT displays a strong downflow appearing just before the formation of the vortex. Note that the displayed 
evolution of the region is longer in time here than in Fig. [T] and the FOV smaller as framed by the black box at 09:04 UT in the same 
figure. 
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Fig. 3. Map of horizontal velocities (FWHM V!0). The length of 
the black arrow at coordinates [0, 0] corresponds to 1.5 km s"^. 
The background image represents the normalized divergence 
field. The solid/dotted boxes extract the FOV covered by the se- 
quences of images in Fig.[TlFig.[2l respectively, as labeled. 



mating the Doppler shift in the line-center- wavelength using the 
center-of-gravity method to the Stokes / profiles (second row). 
For the inversion, only pixels with Stokes Q, U, or V amplitudes 
larger than 3 times their noise levels are taken into considera- 
tion to exclude profiles that cannot be inverted reliably. As ob- 
served in Mgl magnetograms, the two magnetic lobes can be 
clearly distinguished in the \B\ maps displayed in the first row. 
SOT-SP data start earlier than the filtergrams and magnetograms 
described in the previous section. The first 5 frames display the 
evolution of the magnetic structure prior to 08:40, where we can 
see that it is already rotating. The maximum value of the mag- 
netic field strength was close to 600 G at 08:42, though not in- 
cluded in the displayed images. The rms of the errors in \B\ is 
194 G. We defined the errors to be the diff'erences between the 



inferred from the model and the real parameters at optical depth 
log T = -1 . To check the reliability of the inversion, we also com- 
pared these val ues with those obtained by assuming a weak-field 
approximation ("Landi degl'Innocenti 1992) and achieved simi- 
lar results. The longitudinal magnetic flux decreases by about 
30% during the whole time interval. The area of the structure 
seems to shrink with time. In particular, for 09:00-09:20 the 
magnetic patches evolve from being bound together, when the 
two structures lose their identity and become elongated. The Y- 
shape structure that is evident in the Mgl magnetogram (at about 
09:12 UT) is not visible here. One possible explanation is that 
Fe and Mgl lines provide information about difl'erent heights in 
the solar atmosphere. The missing component of "Y" exhibits in 
Mgl magnetograms a very weak and variable signal that is co- 
spatial with bright regions in Call images. However, the compo- 
nent is absent in photospheric filtergrams and SP maps. SP maps 
are available until 09:44 (though not shown in the figure) and 
show that the structure continues rotating after 09:20 when the 
SOT-NFI and BFI time series end, LI and L2 approaching each 
other towards the end of the observing interval. The inclination 
(not shown here) remains stable during the whole observing pe- 
riod. Magnetic flux tubes are vertically oriented with inclination 
values close to 180 degrees. However, changes in the inclina- 
tion can be seen at 08:44, when some BPs appear in-between 
the magnetic lobes. In spite of the high noise level in Stokes Q 
and U, the good agreement between the magnetic structure ob- 
served in the LOS magnetogram and that of SP maps is indica- 
tive of nearly vertical magnetic fields. In the Vlos maps (second 
row of Fig. [2]), the granules are clearly distinguished as upflows 
with velocities from 1 to 2 km s"^ The intergranular lanes are 
the areas showing downflows from -1.3 to -2 km s"^. A strong 
downflow region appears at 08:44 UT (encircled area in Fig.O 
that denotes the location where the photospheric vortex forms 
(Sect. 13.31) . It remains there for several minutes (until 08:56) and 
at 9:20 another downflow region is observed at the location of 
L2 at the end of the sequence, though much less intense. The 
magnetic patches are located in the lanes, with a couple of ex- 
ceptions where the magnetic areas do excursions into granules 
(at 09:00 and 09:07 UT). During the whole sequence, numerous 
small-scale processes take place leading to the formation of BPs 
(Paper I). 
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3.3. Photospheric plasma flows 



The G-band series from SST were used to analyze the hor- 
izontal proper motions of structures in the FOV. Proper mo- 
tions were measured by using the local correlation tracking 
technique (LCT. iNoyember & SimonI 1988) implemented by 
iMolownv^ Horas &Yil (11994 . Maps of horizontal velocities are 
calculated for the time series by using a Gaussian tracking win- 
dow of FWHM r/0, i.e., roughly half of the typical granular size. 
Figure [3] shows the flow map computed from the first of the SST 
time series (si). This map of horizontal velocities is calculated 
by averaging over the total duration of the series (~ 20 min). 
The background image represents in false-color the normal- 
ized (factor of 0.086 s"^) divergence field. Horizontal velocity 
magn itudes are averaging-dependent values ( Vargas Dommguez 
l2008h . though the flow patterns yielded by the analysis of proper 
motions are commonly assumed to represent plasma motions 
thus enabling the detection of general trends, which is our main 
aim here (e.g., sinks or granulation downdrafts where the cold 
plasma returns to the solar interior). The map of horizontal dis- 
placements is dominated by flows coming from granular ex- 
plosive even ts and commonly associated with mesogranulation 
(iRoudier & Muller 2004; Bonet et al. 2005). We identified only 
a pair of examples where strong sinks are present at intergranular 
lanes (coordinates [6.5, 3.5] and [8.5, 7.5] in Fig. [3]) but, in con- 
trast to other areas displaying downflows, these are characterized 
by being the draining point where all horizontal velocity vec- 
tors in the neighborhood converge. This behavior was observed 
in numerical simulations by lNordlundl (1985), who showed that 
the plasma on granular scales in strong downflow regions tend to 
rotate around the center, like a bathtub efl'ect. Both vortex-type 
events last for several minutes (> 20 min). Nevertheless, the one 
framed in the two boxes in Fig.[3]is not only less stable and has 
a defined onset in our data (at -08:48 UT, Paper I) but also ap- 
pears to be the final destination of the magnetic concentrations 
being dragged from an upper location (at coordinates [5, 7] in the 
same figure), as described in Sect. l3.1[ The latter vortex remains 
rather stable but disappears abruptly in the map of horizontal ve- 
locities computed for s2, which is dominated by an organized 
flow coming from the lower right part of the FOV (Fig.|4l right). 
Unfortunately, as already mentioned, there is a 7 -min gap be- 
tween ^7 and s2 meaning that it is not possible to define exactly 
when and how the vortex disappears. 



4. Discussion 

In Fig. 131 we plot the trajectories followed by the cen- 
troid of the distribution of magnetic flux density for lobes 
LI and L2. The centroids denoted with plus-signs and as- 
terisks for the respective lobes are now calculated for the 
individual maps obtained from SP data with a cadence of 
~36 s using the magnetic flux density M for weighting. They 
are defined by Xc = Z^Ii -^/^U/. J/)/ Z ^U/. J/) and yc = 
Yli=iyi^i^i^yi)l Y^^i^hJi)^ where (x/, ji) are the coordinates 
of the /th pixel for each lobe. In spite of its lower resolution, 
we use the SP dataset because its time series extends until 09:44 
UT, allowing us to track the dynamics of LI and L2 over a longer 
period. The trajectories are superimposed on both the respective 
flow maps computed for the SST time series (^7 and s2) and the 
divergence field for each series. The comparison with the cen- 
troid positions calculated using the Mgl magnetograms (top row 
in Fig. [T]) allowed us to estimate the errors to be Exi = ±0.35'', 
€yi = ±0.38'' for LI and e^i = ±0.32", €y2 = ±0.28" for L2. 
From the centroid positions, it was possible to estimate the hor- 
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Fig. 4. Temporal evolution of the magnetic centroid for both 
lobes (Li and L2) in Fig. [T] The trajectories followed by the 
centroids are plotted independently and superimposed on the 
flow maps computed for SST/G-band images in time intervals: 
^7 (left) and s2 (right). The background false-color image is 
the divergence field computed from each respective flow map. 
Alphabet letters correspond to UT times as follows: A-08:40, 
B-08:48, C-08:56, D-09:00, E-09:04, F-09:08, G-09:12, H- 
09:21, 1-09:24, J-09:28, K-09:32, L-09:36, M-09:40,N-09:44. 
The FOV is the same boxed in solid line in Fig. [3] 

izontal velocities for each lobe. While LI motion is confined 
to a radius of approximately 1", L2 travels a much larger dis- 
tance. In the time interval prior to 09:14, its motion seems to be 
strongly influenced by the vortex formed in the coordinates [3.5, 
1.5] in the left panel of Fig.lH increasing its velocity from -1.5 
to 3 km s"^ as it approaches this location. After the vortex has 
disappeared, L2 continues its movement with decreasing veloc- 
ity probably influenced by the continuous evolution of the sur- 
rounding granules towards a new sink denoted by the region of 
negative divergence whose center is located at coordinates [1.0, 
2.5] in Fig. |4l (right). Because of the lack of space here, addi- 
tional analysis of centroids velocities will be deferred to Paper I. 

A possible explanation of the observed magnetic concentra- 
tions being dragged towards the vortex center could be that this 
is the result of the stochastic evolution of granules that allow 
the BPs to approach the vortex influence, increase their veloc- 
ities, and eventually /a// into the vortex (surface scenario). On 
the other hand, long-lasting (> 1 h) vertical magnetic flux tubes 
anchored beneath the surface might be directly influenced by the 
vortex action due to some interaction mechanism taking place 
underneath. Although presenting additional evidence to support 
this scenario is beyond the scope of this paper, this possibility 
should not be ruled out. 
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